Electronic structure beyond the generalized gradient approximation for Ni 2 MnGa 
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The stability of the nonmodulated martensitic phase, the austenitic Fermi surface and the 
phonon dispersion relations for ferromagnetic NLMnGa are studied using density functional the¬ 
ory. Exchange-correlation effects are considered with various degrees of precision, starting from the 
simplest local spin density approximation (LSDA), then adding corrections within the generalized 
gradient approximation (GGA) and finally, including the meta-GGA corrections within the strongly 
constrained and appropriately normed (SCAN). We discuss a simple procedure to reduce a possible 
overestimation of magnetization and underestimation of nesting vector in SCAN by parametrically 
decreasing self-interaction corrections. 


I. INTRODUCTION 

During the last decades, the intermetallic Heusler al¬ 
loys family has attracted enormous interest because a 
wide spectrum of the remarkable properties related to 
elasticity, magnetism, and thermodynamics 1-8 . In par¬ 
ticular, ternary Heusler compounds with generic for¬ 
mula X 2 YZ crystallize in L2i-cubic structure in the 
austenitic phase and can undergo a martenisitic phase 
transition (MT) to lower symmetry structure upon cool¬ 
ing. In general, X and Y are 3d transition metals and 
Z is an sp element of III-V group in the periodic table. 
Usually, the martensitic temperature (Tm) depends on 
the chemical composition, and thus on local interactions, 
plastic deformation, and heat treatment protocols 1,2,5 . 
For most compounds, Tm < Tc , where Tc is the Curie 
temperature. As a result, the magnetic degrees of free¬ 
dom are expected to strongly couple with the lattice, es¬ 
pecially when Tm is close to T^ 1,9-11 . 

In 1996, Ullakko and collaborators studied the mag¬ 
netic control over the shape-memory effect in the ferro¬ 
magnetic Ni 2 MnGa Heusler alloy 12 . These authors were 
able to induce reversible deformations of 0.2% by applica¬ 
tion of magnetic fields below 1 T. In 2002, a larger effect 
of about 10% was obtained in similar compounds with a 
slightly different composition 13 . Today, a great effort is 
deployed to search new Heusler compounds with optimal 
properties for sensors and actuators 2,14,15 . 

Concerning the origin of the MT in Ni 2 MnGa 16 , a the¬ 
oretical explanation based on the Jahn-Teller effect has 
been proposed by Fujii et al. 17 from the calculated den¬ 
sities of states (DOS), which is also consistent both with 
another DFT calculation by Ayuela et al. 18 and photoe¬ 
mission experiments by Opeil et al. 19 . However, the ex¬ 
istence of three-layered premartensetic (3M), martensitic 
five-layered (5M), seven-layered (7M) modulated phases 
in addition to the non-modulated tetragonal phase can¬ 
not be explained by the Jahn-Teller effect only. As sug¬ 


gested by several authors (see e.g. Refs. 20 23 ), the forma¬ 
tion of modulated phases is associated with an anomalous 
transverse acoustic TA 2 mode in the [110] direction at the 
wave vector q = ^ [££0] with £ ~ 0.33 and, consequently, 
a softening of shear modulus C'. Thus, the parent cu¬ 
bic phase becomes unstable because the atomic planes 
can shuffle along the [110] direction and the premarten- 
sitic phase transition occurs at Tp « 260 K 24,25 . In the 
case of MT, the transition temperature is Tm ~ 320 K 
while the softening of TA 2 phonon branch is more pro¬ 
nounced 26 , and it is shifted to £ ~ 0.43. The role of the 
phonon anomaly in destabilizing the austenitic phase can 
be explained by the nesting of the austenitic Fermi sur¬ 
face 27 (FS) and the electron-phonon coupling 26,28 . The 
relationship between the premartensitic lattice soften¬ 
ing with the generalized susceptibility singularities pro¬ 
duced by the FS nesting has been reviewed by Katsnel- 
son et al. 29 . 

For magnetic shape memory compounds, some impor¬ 
tant questions regarding this relationship remain to be 
clarified despite a large amount of work (for instance, see 
Refs. 4,19,26,27,30-33 ) because the lack of information on 
correlation effects present in martensitic ferromagnetic 
phases. Nevertheless, despite correlation effects, the FS 
nesting could stiff play an important role in the insta¬ 
bility of the austenitic phase upon cooling. Velikokhat- 
nyi and Naumov 30 have first explored the FS nesting by 
using density functional theory (DFT) within local spin- 
density approximation (LSDA) and concluded that the 
nesting vector parameter £ extracted from generalized 
susceptibility is approximately 0.42 justifying a MT to 
the 5M phase rather than the premartensitic transition 
to the 3M phase. Lee et al 81 with the same method 
observed a nesting vector in agreement with the experi¬ 
mental phonon anomaly 21,22 by renormalizing the value 
of the magnetic moment with the factor of 70%, which 
corresponds to the magnetization at Tp. 

Haynes et al. 27 conducted a comprehensive study 
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combining positron annihilation experiments with linear 
muffin-tin orbital electronic structure calculations within 
LSDA 34 . These authors found that the peaks of the gen¬ 
eralized susceptibility x(q) appear both in spin-up and 
spin-down channels. Siewert et a/. 32,33 examined the ef¬ 
fect of corrections beyond LSDA on x(q) in the frame¬ 
work of the generalized gradient approximation (GGA). 
Their results confirm that the FS nesting is present not 
only for spin-down channel, but for spin-up as well. They 
also noticed new susceptibility peaks along other direc¬ 
tions produced by GGA. Bungaro et a/. 26 performed both 
phonon dispersion and FS calculations demonstrating 
that GGA corrections to LSDA are beneficial in order to 
improve the agreement between theory and experiment. 

Since correlation effects in N^MnGa appear to play 
a crucial role on total energies, DOS, and magnetic 
moments 35,36 , we also expect major corrections beyond 
LSDA and GGA for describing FS nesting and phonon 
instabilities. A simple step forward is the strongly con¬ 
strained and appropriately normed (SCAN) functional, 
which is the most promising meta-GGA scheme due to 
the number of exact constrains fulfilled by the exchange- 
correlation energy 37 39 . In the present study, a SCAN —V 
parametric method is also introduced to measure the ef¬ 
fect of self-interaction correction (SIC). 

The outline of the paper is as follows. Section II con¬ 
tains the computational methods and calculation details. 
Section III is devoted to the discussion of the results of ex¬ 
change correlation corrections on FS, nesting vector, gen¬ 
eralized susceptibility, and phonon dispersion relations. 
The concluding remarks are presented in Section IV. 


II. CALCULATION DETAILS 

The electronic structure calculations were performed 
using the spin-polarized DFT within the projected aug¬ 
mented wave (PAW) 40,41 method implemented in Vi¬ 
enna Ab-initio Simulation Package (VASP). LSDA and 
GGA parametrized by Perdew-Zunger 42 and Perdew- 
Burke-Ernzerhof 43 , respectively, were used to describe 
correlation energy. Effects beyond GGA are described 
with the meta-GGA corresponding to the SCAN imple¬ 
mentation 39 . A parametric study was performed us¬ 
ing SCAN —U method. The k -points within the Bril- 
louin zone were generated using a uniform Monkhorst- 
Pack 44 mesh of 12x12x12. The cutoff energy was set to 
700 eV. The PAW pseudopotentials were generated with 
the following atomic configurations: Mn(3p 6 3d 6 4s 1 ), 
Ki(3p 6 3d 9 4:S 1 ), and Ga(3d 10 4,s 2 4p 1 ). The calculations 
were converged with the energy accuracy of 10 _6 
eV/atom. Conjugate gradient algorithm was used to 
minimize all the residual forces until the convergence cri¬ 
teria of 0.01 eV/A. The phonon dispersions were calcu¬ 
lated with the PHONON 45,46 package. Additional details 
can be found in Refs. 47,48 . 

For the FS modeling, the 51x51x51 Monkhorst- 
Pack /c-grid was used. The generalized susceptibility was 


calculated for both majority and minority spin channels 
following 26,30,31 


( \ /(e m (k)[l - /(e n (k))] 

w £n(k + q)-£m(k) ’ U 

n,m,k 

where /(£ m ( n )) is the Fermi-Dirac distribution function 
and £ m ( k), £ n (k) are energies corresponding to the m and 
n band at the wave vector k. The peaks of x(q) indicate 
electronic instabilities associated with the FS nesting. 

III. RESULTS AND DISCUSSIONS 

To study the electronic structure, we perform the ge¬ 
ometry optimization of austenitic and martensitic phases. 
We obtain the equilibrium lattice constants and magnetic 
moments for LSDA, GGA, and SCAN. These results are 
presented in SM 48 . 

A. Fermi surface and generalized susceptibility 

Figure 1 shows the band structure in the vicinity 
of Fermi level calculated at equilibrium volume of the 




FIG. 1. (a) Spin-up and (b) spin-down bands crossing the 
Fermi level. The dash, full, and dot lines correspond to the 
LSDA, GGA, and SCAN, respectively. 
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FIG. 2. (left panel) Fermi surfaces of minority spin band of NLMnGa. Blue (orange) color corresponds to 63 (64) band in 
Fig.l, respectively, (right panel) Cross sections of the FS at k z — 0.5(27t/o). Red arrows illustrate the nesting vectors. Dot 
lines delineate the boundaries of the first Brillouin zone. 


austenitic phase. The Fermi level is crossed by three 
spin-up bands and two spin-down bands for LSDA, GGA, 
and SCAN. The FS crossings appear to be robust to¬ 
ward the description of correlation but some differences 
can be observed. The most significant difference is found 
along L-U-W and L-K-X paths for the spin-down bands 
as shown in Fig. 1(b). In this case, band 64 calculated 
within SCAN does not cross the Fermi level along the 
path K-r in contrast to LSDA and GGA indicating that 
some parts of the FS vanish. 


Figures 2 (a-c) show the minority spin FS for the 
austenitic phase within LSDA, GGA, and SCAN, which 
consists of two sheets corresponding to bands 63 and 64. 
The results for LSDA and GGA are in good agreement 
with previous studies 27,33 . As a consequence of the in¬ 
creased magnetic moment in SCAN, a remarkable mod¬ 
ification of the minority FS occurs. In particular, the 
FS sheet associated to band 63 expands, while the sheet 
corresponding to band 64 shrinks. Moreover, the dis¬ 
appearance of a FS piece located at the corner of the 
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FIG. 3. Generalized susceptibility x(q) for NLMnGa. The top, middle, and bottom panels correspond to LSDA, GGA, and 
SCAN, respectively. The first, second, and third columns correspond to minority, majority, and total band contributions, 
respectively. 


Brillouin zone is consistent with Fig. 1(b). 

Figures 2 (d-f) present a visual analysis of the FS cross 
section at k z = 0.5(27t /a) in terms of nesting vectors q x 
and q 2 . These vectors correspond to peaks of x(q) along 
the [110] direction. Our results give |q x | « 0.43 and 0.39 
for LSDA and GGA, respectively, which are in agreement 
with Lee et al 31 and Bungaro et al. 26 . The SCAN value 
is Iq-J ~ 0.20 reflecting significant differences in the FS 
nesting properties. 

A more quantitative analysis of FS nesting proper¬ 
ties visualizes the x(q) 2D-maps for the (110) plane as 
shown in Fig. 3. These spectra are presented for the 
minority, majority, and total contributions for LSDA, 
GGA, and SCAN. In order to highlight important fea¬ 
tures along [110] direction, the profiles of x(q) along this 
direction are also illustrated in SM 48 . According to the 
2D-maps, the results look similar for the minority con¬ 
tributions calculated with LSDA and GGA, while the 


majority contribution presents some additional peaks for 
x(q) obtained with GGA. Interestingly, the total contri¬ 
bution does not differ significantly between LSDA and 
GGA. From the corresponding profiles along [110] direc¬ 
tion shown in SM 48 , it appears that LSDA slightly re¬ 
duces the distance between two susceptibility peaks with 
respect to GGA. The actual values of |q x | are 0.435 and 
0.394 in LSDA and GGA, respectively, while for |q 2 | 
the values are 0.550 in LSDA and 0.596 in GGA. The 
GGA results are in good agreement with the previous 
studies 30,31,33 . SCAN considerable alter the 2D-maps of 
Fig. 3. This modification results in a significant separa¬ 
tion of the two susceptibility peaks and the appearance 
of the additional peaks at q 3 and q 4 originating from 
mixed contribution of spin-up and spin-down. Moreover, 
according to the profile shown in Fig. 3, very different 
nesting vectors |q x | and |q 2 | are found 48 . The values of 
their norms are |q x | = 0.263 and |q 2 | = 0.784. These re- 
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suits show that the correlation effects bring new singulari¬ 
ties in x(q), which can explain a more complex landscape 
of the competing nonmodulated and modulated phases as 
in the case of stripes in cuprate superconductors 49 . 


B. Reduction of SIC 

According to several authors 50,51 , SCAN exaggerates 
the magnetic moment in some 3d transition metals 52,53 . 
The deorbitalization of the SCAN potential has been re¬ 
cently proposed to address this problem 50 . Another sim¬ 
ple way to reduce SIC is to use a SCAN-1/ method where 
U is a parameter measuring the excess of SIC. Figure 4 
shows a parametric study of the SCAN SIC reduction. 
When U = 0, we have full SIC, while a finite U measures 
the suppression of SIC. 

As expected the FS shape becomes more similar to the 
GGA case when a larger amount of U is removed from 
SCAN (see SM 48 ). At the same time, the magnetic mo¬ 
ment decreases and approach the GGA value of 4.17 /ib 
when V ~ 2 eV. Concerning x(q), the reduction of SIC 
leads to the merging of the peaks (q!,q 3 ) and (q 2 ,q 4 ) 
when the critical value U = 1 eV is reached. In this 
case, the 2D-map shown in SM 48 becomes similar to the 
GGA 2D-map. From this parametric study, we deduce 
that SIC associated with SCAN corresponds to an effec¬ 
tive GGA+U scheme with U « 2 eV. This value of U 
is significantly lower than the one proposed by Himme- 
toglu et al 36 (U = 5.97 eV) and the one proposed by 
§asioglu et a/. 54 (U « 4 eV for Mn atom). 



-U, eV 


FIG. 4. Nesting vector parameter (left axis) and total mag¬ 
netic moment (right axis) as a function of the parameter —U 
controlling the self-interaction correction. 



Q[S.S. 0 ] 


FIG. 5. Phonon dispersion in the nonmodulated martensitic 
phase for GGA (red) and SCAN (green). The SCAN disper¬ 
sion presents an imaginary excursion near the T point. 


C. Phonons 

In the case of the austenitic phase, Bungaro et al. 26 
showed that GGA gives imaginary phonon frequencies 
at zero temperature. In SCAN, the appearance of the 
additional peaks in x(q) makes this phase even more un¬ 
stable. This can lead to softening of additional acoustical 
modes. 

Figure 5 shows the phonon dispersion in the case of 
nonmodulated martensitic Llo phase. The GGA phonon 
dispersion yields real frequencies for all phonon wave vec¬ 
tors in agreement with the calculation by Zayak et a/. 55 . 
The SCAN calculation reveals unstable mode near the 
T point as shown in Fig. 5. Therefore, we find that the 
nonmodulated phase is not the ground state as also sug¬ 
gested by Himmetoglu et al 36 . These authors claim that 
the ground state is a modulated martensite. Even within 
GGA, Zeleny et al , 56 indicated that the 40 martensitic 
structure is 2 meV/atom below nonmodulated marten¬ 
site solution. 


IV. CONCLUSION 

Coulomb correlation effects help in stabilizing com¬ 
plex nonmodulated phases needed to understand the 
large magnetic shape memory effect in Ni 2 MnGa 57 . 
Our results indicate that SCAN behaves as an effective 
GGA+U scheme with U parameter of about 2 eV. This 
amount of SIC captured by SCAN modifies the nest¬ 
ing properties of the austenitic FS and thereby creates 
new singularities in the electron susceptibility x(q) jus¬ 
tifying a more complex phase diagram. Consequently, 
the nonmodulated martensitic phase becomes unstable 
within SCAN as demonstrated by the softening of some 
calculated phonon modes. Moreover, the total magnetic 
moment increases significantly as in other 3d transition 
metals such iron 53 . This magnetic enhancement has been 
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considered exaggerated in the literature. Therefore, it is 
possible that SIC corresponding to U = 2 eV should be 
renormalized toward the critical value U = 1 eV, where 
x(q) starts to develop new singularities due to correla¬ 
tion effects. From the present results, we conclude that 
the amount of SIC in N^MnGa can be measured by an 


effective U in the interval of 1 — 2 eV. 
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